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electrons in liquid helium, we demonstrate that the decay of Rabi oscillations 
observed in experiments on 9 Be + can be attributed to the polarization phe- 
nomena investigated here. The measured sensitivity of the damping of Rabi 
oscillations with respect to the vibrational quantum number of a trapped ion 
is also predicted in our polarization model. 
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I. INTRODUCTION 

The experimental techniques developed over the last decade for manipulating trapped 
ions, Rydberg atoms in cavity QED, and traveling fields have enabled a fruitful dialog to take 
place between theoretical and experimental physics, resulting in a mastery of fundamental 
quantum phenomena at a level that seems to herald a new phase in the technology of 
communication |T| and computation 0,0]. The improvement of techniques for generating 
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and detecting nonclassical electronic and vibrational states of trapped ions and both trapped 
and traveling nonclassical states of the radiation field have been of growing interest in a 
variety of experimental applications ranging from quantum measurement concepts |4]|| to 
teleportation || and quantum logic operations 0. In its turn, theoretical physics has 
exploited the possibility of engineering nonclassical states in some exotic applications such 
as detection of gravitational waves |§, or for measuring particular properties of the radiation 
field, such as its phase || or its Q function flO| . Laser manipulation of a string of ions in a 
linear trap has been proposed as a way of implementing quantum gates with cold trapped 
ions ||. Furthermore, interfering laser beams have been used to induce arrays of microscopic 
potentials, the optical lattices ||11|| , by the ac Stark effect, opening the way to simulated spin- 



spin interactions between trapped atoms, similar to those characteristic of ferromagnetism 
and antiferromagnetism in condensed matter physics ||12|| . Quantum computation with ions 
in thermal motion has also been suggested |]T3| , p!4| . 

The practical realization of such interesting proposals in quantum communication and 
computation has come up against the decoherence of quantum states, owing to the inevitable 
action of the surrounding environment [T5| and the intrinsic fluctuations in the interaction 



parameters required for logic operations |TB[. In fact, the need for huge superpositions 



of qubit states for such operations imposes the requirements that the quantum systems be 
totally isolated from the environment and that the interaction parameters involved be tightly 
controlled. Hence, investigation of the noise sources in such promising quantum systems 
turns out to be a crucial step toward the implementation of a quantum logic processor. 

Unlike processes in cavity QED, where decoherence is caused by the well-known cavity 
damping mechanisms, spontaneous atomic emission, and inefficiency of the ionization de- 
tectors, in the domain of trapped ions a large range of error sources have been identified. 
The effects of spontaneous emission have been studied fl7|-|T9|1, as well as dephasing due to 



the ions' zero-point motion [50]. Alternative sources of decoherence have been introduced 



phenomenologically []2T| and stochastic models have been proposed to handle intensity and 
phase fluctuations in the exciting laser pulses, in addition to fluctuations in the ion trap 



potential f22|| . Recently, instead of a stochastic mechanism, Di Fidio and Vogel J2ij have 



put forward a model in which the observed damping in Rabi oscillations |Mj is caused by 
quantum jumps to an auxiliary electronic level. Other significant sources of error are the 
instabilities of the trap drive frequency and voltage amplitude P|25|j. 

Collisions of the trapped ion with residual background gas (usually H 2 in the NIST 
9 Be + experiments |[26|| ) can also be an important error source, although experiments are 
typically carried out in a high-vacuum environment |26| at a pressure around 10 -8 Pa. In 
connection with this, the current article discusses the polarization effects on the residual 
background gas induced by the oscillating trapped ion. We demonstrate that, apart from 
other known processes, the damping and heating mechanism of trapped ions is produced by 
local polarization of the residual thermal background gas (BG). Since the density of the BG 
is rather low, around 10 6 cm -3 p6 |, we do not expect to find a quasiparticle consisting of 
the trapped ion together with its surrounding polarization cloud, the so called polaron. The 
density of the BG does not permit a polaron binding energy and the trapped ion is scattered 
by the BG. For practical purposes, in our model we assume that the BG is continuous and 
that the ion is scattered by its surface oscillations, by analogy with the interaction of surface 
electrons in liquid helium (where the oscillations are called ripplons in their quantized form) . 

In contrast to loss mechanisms revealed by damped Rabi oscillations, it is well known 
that background gas can heat trapped ions by transferring energy during an elastic collision. 



A heating rate can easily be estimated from the total collision cross section [pqj . Although 
elastic collisons are expected to be the main source of decoherence in ionic traps, inelastic 
collisons also take place, changing the internal state or even the species of the trapped ion. 
In experiments on 9 Be + , inelastic processes may convert the ion to BeH + (upon collision 
with a H 2 molecule) when resonant light is applied to the 2 Si/ 2 — > 2 Pi/2,3/2 transitions 0. 
Both types of inelastic collision, chemical reactions and charge exchange, besides depending 
critically on the constituents of the BG, occur only when the interparticle spacing between 
the trapped ion and the neutral background atoms approaches atomic dimensions. In our 
model, a mean interparticle spacing between the ion and the BG atoms is introduced. An 



upper limit on the rate of inelastic collisions follows from the Langevin rate, which accounts 
for the background neutrals that penetrate the angular momentum barrier and undergo a 
spiraling-type collision into the ion p6| , p7|j . 

The attractive interaction potential resulting from polarization of the neutral background 
by the electric field of the trapped ion is given by U(r) = ~XQ 2 / (87reor 4 ) , where \ is the 
polarizability and q the ionic charge. Spiraling collisions result when the impact parameter 
is less than a critical value p = (x<? 2 / 'vreotriD 2 ) 1 ^ 4 : , where m and t> are the reduced mass and 
relative velocity of the pair. From the critical impact parameter follows the Langevin rate 
constant £ = 7rp 2 D, which gives to the reaction rate 1Z = pC = pq (rr x/eom) 1 ^ 2 , where p is 
the density of the BG. From the parameters used in the NIST 9 Be + experiments, we obtain 
a small estimated probability of inelastic collisions with the BG constituents. 

In addition to the chemical reactions and charge exchange, the BG can heat or cool the 
trapped ion through energy transfer during an elastic collision. An estimate of the elastic 
collision rate can be inferred from the total collision cross section a in a A/r 4 potential. 
Assuming A = U(r), it follows, from a conservative estimate, that elastic collisions will 
also be rare [^(J. However, the effects of collisional heating (cooling), as well as the model 
we present here, can be tested, as suggested in Ref. [ Z6fl , by raising (lowering) the BG 



pressure. On the other hand, it is well known that when ions are first loaded into a trap, 
elastic collisions with the BG are beneficial, allowing laser cooling to proceed faster, the 
BG providing a viscous medium and bringing the temperature of the trapped ions into 
thermal equilibrium with the surrounding atoms. We conclude from the work reported here 
that the damping medium provided by the background has also to be taken into account 
when the trapped ions are cooled to their motional ground state. This claim is supported 
by the excellent fit between our model and the experimental data for measurement of the 



fluorescence probability of the electronic ground state, reported in ||24 . 

Finally, we mention that, for a linear Paul rf trap, the process of decoherence may be 
dominated by that of the motional state, instead of those due to internal levels, caused by a 



nonideal applied field ||26|| . In fact, the internal levels of the trapped ion are metastable, typ- 



ically two ground state hyperfirie sublevels f28fl . While Refs. |2^] take into account technical 



fluctuations in the applied fields used to manipulate the trapped ions, Ref. [^3|] attributes 
the error source to the coupling of the internal states to the environment. In this paper we 
present a model that accounts for the dominant decoherence of the ion motion in terms of 
the induced polarization of the BG. 

The paper is organized as follows. In Sec. II we present a detailed discussion of the 
fundamental ion-laser and ion-BG interactions. In Sec. Ill we examine the effects of the 
ion-BG coupling on the behavior of a prepared motional-electronic state of the trapped 
ion, when considering the Carrier and the anti-Jaynes-Cummings ion-laser interactions. For 
the latter ion-laser interactions the Rabi oscillations of a trapped ion initially cooled to 
its motional ground state are computed under the ion-BG coupling, and our results are 



compared with the available experimental data for a trapped 9 Be + ion p4j. Sec. IV is 



dedicated to comments and conclusions and, finally, in the Appendix A we show how to 
obtain a Frohlich-type ion-BG interaction, following the model used for surface electrons in 
liquid helium. 

II. ION-LASER AND ION-BG INTERACTIONS 

We consider a single trapped ion of mass m in a one-dimensional harmonic trap of fre- 
quency v. The ion has forbidden transitions between two internal electronic states (excited 
| X) an d ground | J.) states, assumed to be hyperfine sublevels of the ground state), separated 
by frequency ujq and indirectly coupled by the interaction with two laser beams, with fre- 
quencies u>i and w 2 , in a stimulated Raman-type configuration. As indicated in Fig. 1, the 
laser beams are detuned by A from a third more excited level |r) which, in the stimulated 
Raman-type configuration, is adiabatically eliminated when A is much larger than three 
quantities: the linewidth of level |r), the coupling associated with the | f) ^ l r ) and | |) 



|r) transitions, and the detuning 5 = luq — lul (lul = uj\ — UJ2) p6|j28|j29|] ■ The ion inter 



acts with an effective laser plane wave propagating along the x direction, with wave vector 
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Ivl = ojl/c. In this configuration, only the ionic motion along the x axis will be modified. 
The transition between | j) and a fourth level \d), achieved by another laser strongly coupled 
to the electronic ground state, is analyzed in order to measure the ionic vibrational state by 
collecting the resonance fluorescence signal, which is the probability of the ion being found 
in the internal state ||) J3"D| . 



The ion-laser interaction Hamiltonian that describes the effective interaction of the quan- 
tized motion of the ionic center of mass (CM) coupled to its electronic degrees of freedom 

is 

H lon . laser = Ml (<t + e^—^ +a_ e -^+^-^) , (1) 

where a + = | f)(J, |, cr_ = | | and a z are the usual Pauli pseudo spin operators, x is the 
position operator for the x coordinate of the ion, Q is the effective Rabi frequency of the 
transition ||) <-> ||), and is the phase difference between the two lasers. 

The ion-BG interaction (the polarization of the BG induced by the oscillation of the 



trapped ion) will be described by a Frohlich-type electron-phonon Hamiltonian ||31|| . In 



Appendix A we show, by analogy with surface electrons on liquid helium []52|-p5|], that a 
polaron-like interaction results when the electric field of the trapped ion polarizes the neutral 
BG (usually H2 in the NIST 9 Be + experiments f26|). The attractive ion-BG interaction 
potential is given by U(r) = —XQ 2 / (8vreor 4 ), where x is the polarizability and q the ion 
charge, and the resulting Hamiltonian reads 

Hion—BG = Wk { bk ^ +b l ' ( 2 ) 

k 

where b\ (6^) is the creation (annihilation) operator of BG-oscillation quanta, Vk stands for 
the coupling strength, and k are the x-components of the BG-oscillation wave vector. In a 
frame rotating at the "effective laser frequency" ool, the ion-laser and ion-BG Hamiltonians 
are given in the Schrodinger picture by (h = 1 from here on): 

H lon . laser = n (a+ e^ a+a ^ +a. e -^ a+a ^ +i ^ , (3) 
Hion—BG = J2 Vk { hk e l%(aW) +b{ e"^( a+at )) , (4) 



where at (a) is the creation (annihilation) operator of vibrational quanta, tjl = k^j \/1mv is 



the Lamb-Dicke parameter, and r\ k = k/yZrrw stand for Lamb-Dicke-like parameters due 
to the ion-BG interaction. The total Hamiltonian is given in the Schrodinger picture by 
H = H + H ion _i aser + H ion _BG, where H indicates the free Hamiltonian composed of the 
internal and motional degrees of freedom of the trapped ion plus the BG: 

x 

H = va ] a + -a z + 2J u k b ] k bk- (5) 

k 

Writing the total Hamiltonian in the interaction picture (bold labels), by the unitary trans- 
formation U(t) = exp (—iH t), and then expanding the resulting expressions in terms of the 
parameters rji and r] k , we get 

\m,l=0 ' ' / 

Kion-BG = £ V*** ( £ {jlI ^ha^a l ^-i)^ h] t +hc \ (7) 

k \m,l=0 ' ' J 

where ion-laser resonance is achieved by tuning the laser frequencies to obtain 5 = —Iv 
(£ = m — l). In what follows some reasonable approximations are made, in order to simplify 
considerably the Hamiltonians ([|) and (0). First we mention (z) the standard Lamb-Dicke 
limit for the ion-laser interaction, for which t]l <C 1 (where the ionic CM motion is strongly 
localized with respect to the laser wavelengths). Next, from the low energy of the BG 
oscillations (ii) a Lamb-Dicke-like limit, % < 1, will be assumed for the ion-BG interaction. 



This limit is somewhat analogous to the case of large polarons in solid state physics [|36 
In fact, the motion of large polarons is continuous, as should be that of a trapped ion. In 
contrast, small polarons recognize the periodicity of a solid, becoming localized as in the 
strong coupling theory and assuming atomic dimensions. With these two approximations 
we obtain the simplified Hamiltonians 

H lon _ laser = n (a+e-^+z^a+ate^^'-^+z^a+ae-^^^^+^.c.) , (8) 
Hion-BG = J2 Vk ( bk e ~ iUkt + ir lkbk^ e^-"^ +i Vk b k a e-^+^ +h.c.) . (9) 



With the addition of (Hi) the optical rotating wave approximation, three specific Hamil- 
tonians for the ion-laser interaction are obtained, depending on the choice of the ion-laser 
detuning: 

(a) the carrier Hamiltonian (5 = 0), 

HL- ittser = fiKe-na_e^), (10) 

which induces the transition \n, I) < — ► \n, t) (where \n) indicates a motional Fock state), 
and is responsible for rotating only the internal electronic levels of the ion wave function in 
accordance with 

e - tH ion-laser T | n? — COS (Qt) \n, |) ~ 1^ Sm (Qt) \n, |) , (lla) 

e -iHg n _ laser r | n> ^ = cog ^ | Wj ^ _ ie -iv gin ^ | Wj ^ . ( Ub ) 

(b) the Jaynes-Cummings Hamiltonian (8 — v), corresponding to the first red sideband, 

^on-iaser = *fcfi {a+ae^ - a^e +i *) , (12) 

induces the transition \n, j) < — ► \n — 1, ]), in such a way that the electronic and vibrational 
modes evolve as 

e -^f c n _ laser r | n? T) = ^ | n> T) _ e -i? Sn | n + 1? ^ ) (13a) 

e -^-«..er^ | n , |) = C n _! |n, I) + e^S n ^ \n - 1, T> , (13b) 

where C„ = cos(gry/n + 1), = sin (pr\/n + 1), r is the duration of the laser pulses, and 
9 = Vl^; 

(c) the Anti- Jaynes-Cummings Hamiltonian (5 = —v), corresponding to the first blue 
sideband, 

H£t C -iaser = ^ [a + a^ - a _ae^\ , (14) 

induces the transition \n, I) < — > \n + 1, |), and the electronic and vibrational modes evolve 
as 
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e -u^ w r | n> |) = Cn | nj |) + e -^ Sn | n + 1,1), (15a) 
e-^w^ | n , |) = C n _x |n, T) - e^SU \n -1,1). (15b) 

Finally, assuming hypothesis (Hi) and that (iv) the BG-oscillation modes are closely spaced 
in frequency, with the maximum of their spectrum far from zero, the ion-BG Hamiltonian 
becomes 

Hion-BG ^i^VkVk [ha^e^-^ -ft+ae-^-"*)) , (16) 
k 

which plays the role of a reservoir for the motional degree of freedom of the trapped ion. 

Next, we analyze the effects of the ion-BG coupling on the time evolution of the ionic 
external and internal degrees of freedom. We note that the ion-laser detunings 5 = ±z/ cause 
the internal levels of the trapped ion to be affected by the BG, through the motional degree 
of freedom, differently from the choice 5 = 0, which leaves the BG without effect on the 
internal levels of the trapped ion. 

III. EFFECTS OF THE ION-BG COUPLING ON A PREPARED 
MOTIONAL-ELECTRONIC STATE 

In what follows, the effects of the ion-BG coupling on the behavior of a prepared motional- 
electronic state of the trapped ion will be considered for the ion-laser detunings 5 = and 
—u, representing the Carrier and the Anti-Jaynes-Cummings ion-laser interactions discussed 
above. 

A. Carrier Hamiltonian 

First, using Glauber's P representation [K7J] we briefly analyze the dynamics of the ionic 
motional state when 5 = 0, the total Hamiltonian of the system, in the Schrodinger picture, 
being 
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H = H + U\t)HZ n _ laser U(t) + U\t)H lon . BG U(t) 



vera 



+ "Ah + n ((T+e-* + <T-e**) + % ^ rj k V k (b k a? - b\a>j (17) 



Focusing on the time evolution of the ionic motional state, which is not affected by the 
internal levels in the Carrier regime, we observe that the Hamiltonian ( |TTD is analogous to 
that of a single-mode field trapped in a lossy cavity. The latter Hamiltonian has received 
considerable attention recently for computing the fidelity of a cavity-field state required 
for quantum communication or computation purposes In this connection, we find it 



interesting to analyze the decoherence process of a prepared motional state of a trapped ion 
due to its coupling to the BG. 

Solving the coupled linear equations of motion for the operators a and b k , we obtain from 
Eq. QT7D 

a(t) = u(t)a(0) + v k (t)b k (0), (18) 
k 

where, disregarding the typical small frequency shifts and introducing the damping constant 
T, the time-dependent coefficients, obeying the initial conditions u(0) = 1 and v k (0) = 0, 
can be written as 

u (t) = e -( r / 2 +-)', (19a) 

1 _ e ik-")f e _r */ 2 
v h (t) = -VkV k -. (19b) 

r/2 - i(u k - v) 

From the above results, we obtain the normally ordered characteristic function Xat(£>£)> 
defined in the Schrodinger picture as the expectation value 

XN fa t) = Tr \p(t) e« at W e ra Wj , (20) 

where the density operator p(t) refers to the composite ion-BG system, assumed, as usual, 
to be initially decoupled, p(0) = pionifypBaiQ) (the ion-BG coupling is turned on suddenly 
at t — + ). For a motional state of the trapped ion, prepared as Yli c i \ a i)-> follows the 
representative term of the reduced density operator /w(0) = {oli\, while for an initially 
thermal state with Gaussian distribution the reduced density operator for the BG reads 

10 



_ r e -\M 2 /{n k ) 

(°) = ny n{nk) iMWh. (2i) 



Pbg 

k 

From the characteristic function f^O!) we derive the conditional distribution function 



P (^) = ^ J ^-^ XN (Z,t)d 2 Z 

(7*-«*(t)at)(7-«(t)a 2 ) 



(ai|a 2 ) 



w) J ■ (22) 

whose dispersion is given by D(t) = J2k ( n k) l^feC^) | 2 = ^th (l — e_r< )' where the occupation 
number of the BG-oscillations in thermal equilibrium reads n t h — (e _l///fcsT — l) , ks is the 
Boltzmann constant, and T is the absolute temperature. Fig. 2 display the damping process 
in the evolution of the conditional distribution function P(7, t) for a prepared motional 
coherent state. The rotation in phase space arises from the factor exp(— ivt) in Eq. ( |19a| ) and 
we observe both the effects coming from the environment: the loss of excitation, carrying the 
initial coherent state to the vacuum state, and the diffusion due to the nonzero temperature 
of the BG. These effects are clearly seen from the shape of the distribution function P(j,t) 
depicted in Fig. 2 at Tt = 0.2 and 0.9. 

From the distribution function P(j,t), we can obtain the reduced density operator for 
the ion, 

Pion{t) = ^ J P{j,t) |t) ( 7 | A, (23) 

and the mean value of operators associated with the ionic motional states; for example, the 
mean excitation of a prepared motional vacuum state, i.e., ot\ = a 2 = 0, reads 

<a f (t)a(t)) = J P(a,t)a*ad 2 a = D(t), (24) 

showing the heating process of the motional vacuum state due to a thermal BG. 



B. Anti-Jaynes-Cummings Hamiltonian 



Next, we examine the action of the ion-BG coupling on a prepared motional-electronic 
state of the trapped ion when 5 = —u, the total Hamiltonian being 
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H = H + U\t)Hi J n c _ laser U(t) + tf(t)H ion - BG U(t) 

= i>a)a - ^er z + ^ uJkb\b k + ir]L& (a+a)e~ lv - cr^ae^) + i ^ r] k V k (b k a) - b\a^ . (25) 

k k 

We are considering the Anti-Jaynes-Cummings (instead of the Jaynes-Cummings) Hamilto- 
nian for the ion-laser interaction, so as to compute, in this regime, the damping of the Rabi 
oscillations of a trapped ion initially cooled to its motional ground state and to compare our 
results with the available experimental data for a trapped 9 Be + ion |24 . 



To compute the atomic inversion from Hamiltonian (^), we use the techniques de- 
veloped to investigate the Jaynes-Cummings Hamiltonian for an atom interacting with a 
monochromatic field trapped in a lossy cavity at a finite temperature, which is analogous 
to Eq. ([25|) and has received considerable attention in the literature |38|- f41|1 . The standard 
master-equation technique was employed in Ref. |39| to obtain the density matrix for the 
combined atom-field system in the interaction picture. However, following [J38| , we calculate 
the evolution of the atomic inversion from the time-dependent averages of atomic excitation, 
photon numbers, and other field quantities. Although also based on the master-equation 



approach, the technique in consists of working with a hierarchy of c-number quantities 
derived from the equation of motion for any operator of the form (a^)™ cl u Oa (where Oa is 
an atomic operator), given by 

j t [(<f) m a n O A ] = -i [(a^ m a n O A ,~a z +i VL n(a + a^e-^-a.ae + ^) 

where the last term includes the commutation relation between (a^) m a n OA and the compo- 
nents of the Hamiltonian (^) describing the energy of the BG and its interaction with the 
ionic external modes, as well as the motional energy of the trapped ion. Under the Wigner- 
Weisskopf (WW) approximation and supposing that the BG is in thermal equilibrium, we 
obtain 

(Jfc [faT "]) = (iv{m-n)-~{m + n?j ((^) m a n ) BG 

+Tmnn th ^) m ' 1 a n ' 1 ^ BG . (27) 
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The damping constant, arising from the WW approximation, is T = 2tt [g(z/)] 2 A(z/), where 
g(u) = g v / c is the coupling constant evaluated at k = vjc and A(u) = W^/ttc 3 (V being the 
quantization volume) is the density of states. The occupation number n t h follows from the 
assumption of thermal equilibrium, where the noise operators M{t) = Ylk 9 k ^k{^) e ~ l ^ k t 
(with g k = r) k V k ) satisfy (N{t)) BG = {U\t)) BQ = 0, and (jVt(t)jV(t')) BG = Tn th 5{t - t') 
(see Ref. @). 

From Eqs. (|26D and (|27|) we derive the time-evolution of the atomic inversion (a z ) and 
the motional occupation number of the trapped ion : 

d (<7 Z ) 



dt 
d (a^a) 



2g (a+ate-* + a.ae +ilp ) , (28) 
g (a+aV^ + a.ae +ilp ) - V {a) a) + Tn th (29) 



dt 

with the angular brackets denoting the BG as well as the quantum mechanical average. 
Eqs. (EH) and (|29|) involve the average of the Hermitian operator a + a^e~ tip + a_ae +llp whose 



equation of motion involves the quantity (a* adz}. In its turn, the equation of motion 
of this quantity involves higher-order products of averaged operators such as ^(a^) 2 a 2 ^, 
^(a^) 2 a 2 cr 2 ^, and so on. To deal with such averaged operators it is convenient to define the 
c-numbers variables 

V n =((a?) n a n ), n>0, V = 1; (30a) 
Qn =((<#)" a n a g ), n>0, Q = (a z ) ; (30b) 
TZ n = (a + {a)) n a^e^ + a_{a)) n ~ l a n e +iip Y n > 0. (30c) 

Solving the dynamical equations for the above c-number variables we obtain the dynamics 
of the atomic inversion and the motional occupation number of the trapped ion. From 



definitions (|30a|) - (|30c|) , and Eqs. fl26|) and (|27|) we obtain the following set of equations 
dV 

= nglln - nYV n + n 2 Yn th V n _ x , (31a) 

dO _ 

ngll n + 2gK n +i ~ nTQ n + n YrhhQn-u (31b) 



dt 
dTZ n 



2gQ n + ngVn-i - ngQn-i - (n - l/2)m n + n{n - l)Tn t hR-n-i- (31c) 



dt 
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Note that at zero temperature (n t h = 0) and an initially prepared ionic state \n = 0, J,) 
we obtain a closed set of equations, since the expectation value of the operators involving 
quadratic or higher powers of the ionic motional operators a and a) are therefore zero 
at all times |[41]| . However, to compare the decay of Rabi oscillations resulting from our 
model with the experimental data [|24[, we have to deal with the realistic case of non-zero 



temperature. Observing that for the initially prepared ionic state |n = 0, j), the Anti- 
Jaynes-Cummings Hamiltonian induces the transition |0, \) < — > |1, |), and it is possible to 
estimate a temperature around v/ks for the trapped ion and, consequently, for the BG. In 
fact, under such conditions the thermal energy of the trapped ion fc^T is around z/, resulting 
in a motional occupation number around unity (n t h ~ 1). On the other hand, since the 



Rabi oscillations in NIST experiments [24 survive considerably above 100/is, it is possible 
to estimate a damping factor V around 10 _3 s. Having these values at hand and truncating 



the system of equations ( |31a|) - (|31c|) at n = 4 (an approximation due to the low motional 



occupation number estimated above for the BG), we obtain, via the numerical Laplace 
transform method, the time-evolution of the probability for fluorescence measurement of the 
electronic ground state Pj.(t) = (1 — (<J Z )) /2. Fig. 3 displays the behavior of the function 
Pl(t) computed from our model, which is in excellent agreement with the experimental 
data reported in ||24|| . We used in our numerical calculation the experimental parameters 
t]l = 0.202 and Q/2n « 475 kHz. In particular, our model reproduces the asymmetry of 
the decay of Rabi oscillations seen in the experimental data. To obtain such behavior, we 
used the parameters n t h = 1.0 and T/g = 6.0 x 10~ 3 , which are in good agreement with the 
values estimated above. It is worth noting that, when switching off the laser pulse, Eq. (|24"D 
indicates that the motional occupation number of the trapped ion {Va) goes asymptotically 
to unity, as expected. 

In order to demonstrate from our model the sensitivity of the damping of Rabi oscillations 
to the motional quantum number, we depicted in Figs. 4(a,b) the behavior of the function 
P^(i) computed from our model (full lines) when the initially prepared ionic states are 
\ n — 0, |) [Fig. 4(a)] and \n = 1, 1) [Fig. 4(b)]. Figs. 4(a,b) also display the heuristic 

14 



relation used in |24]] for fitting the experimental data, 

' (32) 



where p n is the initial probability distribution of the motional states in the Fock basis and 
In = 7o{n + 1) ' 7 is a phenomenological damping rate. The behavior of the heuristic relation 



(|32|) for | ra = 0, |) (Fig. 4(a)) and \n = 1, j) (Fig. 4(b)) is shown in dotted lines using the 
value 70 = 11.9(4) kHz estimated in |[24| . We stress that Eq. (|32l) describes a symmetric 
decay of the Rabi oscillations instead of the experimentally observed asymmetry which is 
in agreement with our model (see Fig. 3). Therefore, we claim that the polarization of the 
residual background gas is the mechanism leading to the n-dependence of the decoherence 
effects of a trapped ion. 

IV. COMMENTS AND CONCLUSIONS 

We have proposed an alternative mechanism to explain the damping of Rabi oscillations 
of a 9 Be + ion initially cooled to its motional ground state. Our model allows for the influence 
of the polarization of the residual background gas induced by the oscillating trapped ion. 
The polarization of the BG, which in turn influences the trapped ion motion, is described by 



a Frohlich-type ion-BG interaction with a Lamb-Dicke-like limit where % = k/y2rrw <C 1, k 
being the BG oscillation wave vectors, and m (v) being the mass (frequency) of the trapped 
ion. 

Since experiments on trapped ions are carried out in a high-vacuum environment, the 
density of the BG, around 10 6 cm -3 | 2Efl , does not permit a polaron binding energy and 



the trapped ion is scattered by the BG. Both elastic and inelastic collisions are present, 
although elastic collisions are expected to be the main source of decoherence. The trouble- 
some inelastic processes of chemical reactions and charge exchange are expected to be rare. 
Although a conservative estimate for the elastic collisions rate also leads to a small value, 
we claim that the damping medium provided by the background, which allows laser-cooling 
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to proceed faster when ions are first loaded into a trap, has also to be taken into account 



when the trapped ions are cooled to their motional ground state. As mentioned in ||26|| , for 
a linear Paul-rf trap the main source of decoherence may be associated with errors arising 
from the external degrees of freedom of the trapped ion rather than internal levels or a 
nonideal applied field. 

Under the Lamb-Dicke-like limit adopted here, the resulting ion-BG interaction consist- 
ing of the standard "rotating wave" terms (assuming that the spectrum of the phonon modes 
presents its maximum far from zero) turns out to be the usual linear response model of the 
reservoir (the BG) to the system (the trapped ion) |37|,f42|. Glauber's P representation was 



used to analyze the dynamics of the ionic motional state when considering a Carrier pulse 
for the ion-laser interaction. The Carrier pulse prevents the BG having any effect on the 
internal levels of the trapped ion and either the damping or the heating process of the ionic 
motional state, depending on the temperature of the BG, was investigated. On the other 
hand, a master-equation approach was used to achieve a theoretical fit to the exponentially 
decaying, sinusoidal Rabi oscillations [f23j] resulting from a prepared motional-electronic state 
of the trapped ion under the anti-Jaynes-Cummings pulse. 

The complete agreement between the measured data for Rabi oscillations with the be- 
havior computed from our model indicates that polarization effects constitute a relevant 
source of error in the ionic trap. Evidently, the error from polarization of the BG has to 
be taken together with the errors coming from the internal levels and the application of 
a nonideal field p6| , as investigated in Refs. Like the stochastic approach in fl23fl . 



our model also reproduces the asymmetry of the decay of Rabi oscillations, as observed 
in the experimental data. Since the authors in Ref. |S| claim that the asymmetry of the 



measured data is consistent with the assumption of the mechanism of quantum jumps, we 
mention that it is also consistent with the damping-heating competition produced by the 
ion-BG interaction (through the motional occupation number n t h and the damping factor 
T). Moreover, the consistency of our model is confirmed by the fact that the best fit of the 
measured data was accomplished with parameter values (n t h = 1-0 and T/g = 6.0 x 10 -3 ) in 
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agreement with theoretical estimates. 

Finally, we have demonstrated with our model that the damping of Rabi oscillations 
is sensitive to the motional quantum number, as experimentally observed. The behaviour 
of the curves in Figs. 4(a,b), representing two ionic states initially prepared with different 
motional quantum numbers |n = 0, |) and |n = 1, |), respectively, strongly suggests that the 
polarization of the residual background gas is the mechanism leading to the n-dependence 
of the decoherence effects of a trapped ion. 
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Appendix A 

In this appendix, reasoning by analogy with the physics of surface electrons on liquid 
helium [|53| , |S5[| , we show how to obtain a Frohlich-type ion-BG interaction. For practical 
reason, we assume a continuous BG and a cylindrical symmetry for the ion-BG system, as 
depicted in Fig. 5. The ion, oscillating along the x-axis, is assumed to be at a distance z 
from the origin of the coordinate system which is fixed in the center of the BG surface located 
in the xy-plane. This distance z (which contributes to the coupling parameter V&) is here 
assumed to be a mean distance between the ion and the atoms composing the BG, a mean 
interparticle spacing between two colliding partners. As mentioned in the Introduction, the 
ion is scattered by the oscillations of the BG surface, whose dynamic roughness is given by 
£ (P'^j (see Fig. 5). 

The potential energy of the ion due to polarization of the BG is 

XQ 2 A?) . , 1 



u(r,z) = / dV / v ' dz' 







2 
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Substituting the variables z — £ (j^J — C> z ' ~ £ (^J = C> an d defining A/n = %pe 2 /2, we 
can rewrite the potential energy as 



U(~f,z) =5U(T,z)+U Q (r,z) 



5U (T, z) 



U (T,z) 



+ (C - C) 2 



A 

7T 



2 ' 



1 



-i 2 • 



+ (C - C) 2 



Assuming that the roughness of the BG, £ f^) — £ * s sufficiently small and expanding 
£ (f^J an d £ ("r*) in Fourier series we obtain 



5U(T,z) 



7T 



A 



+ (z' - zf 

d 2 t 



z'=0 



k 



(X 2 + z 2 Y J (X 2 + z 



d 2 te ikX 



where S is the area of the BG surface and we have substituted the variable r — "r* = A . 
After integrating over A we obtain the perturbed potential energy 



k -~r* 



- fcKi (Jfez) 



where (kz) is the modified Bessel function. At this point we can quantize the dynamic 
variable £ ( k ^ as follows [!J3|: 



C- 



\ 2p, 



w ( k 



where u [ k ) is the dispersion relation. Defining the coupling parameter 



AC 



t r 1 



V- 



k 



Vs [z 



JfcKi(fcz) , 



z 



we finally obtain the ionic quantized potential energy due to interaction with the BG as 



which results, in one dimension, in the interaction Hamiltonian (|2]). Uq (It*, z) plays the role 
of a reference energy. 
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Figure Captions 

FIG. 1. Electronic energy level diagram of a trapped ion interacting with laser beams 
of frequency uj\ and co> 2 , where 5 = u — uj\ + uj 2 (5 <C A), |r) (adiabatically eliminated) 
is an auxiliary electronic level which indirectly couples the levels ||) and ||), and \d) is an 
electronic level used to measure the fluorescence emission. 

FIG. 2. Conditional distribution function P(^y,t) for an initially prepared motional 
coherent state depicted at Ft = 0.2 and 0.9. The mean value of the amplitude moves on an 
exponential spiral displaying both the effects due to coupling with the environment: the loss 
of excitation, which carries the initial coherent state to the vacuum state, and the diffusion 
due to nonzero temperature of the BG. 

FIG. 3. Time evolution of the probability of fluorescence measurement of the electronic 
ground state Pj.(t), computed from our model (full line) and measured in the NIST 9 Be + 
experiments |2^] (dots). The parameters used in the experiments rjL = 0.202 and Q/2n « 475 
kHz, were adopted for the numerical calculation, together with the values nth = 1-0 and 
T/g = 6.0 x 10~ 3 which are in good agreement with theoretical estimates. 

FIG. 4. Time evolution of the probability of fluorescence measurement of the electronic 
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ground state Pj,(i) for the initially prepared ionic states (a) |n = 0, j) and (b) |n = 1, |), ac- 
cording to our model (full lines) and the heuristic relation (|32| ) (dotted lines), demonstrating 
the measured sensitivity of the damping of Rabi oscillations to the motional quantum num- 
ber 0. 

FIG. 5. Sketch of ion-BG interaction, showing the ion oscillating in the x-direction, 
assumed to be at a distance z from the origin of the coordinate system which is fixed in the 
center of the BG surface located in the xy-plane. The ion is scattered by the oscillations of 
the BG surface, whose dynamic roughness is given by £ (r ) • 
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